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NOTATION 


a o 

B(X,t,T*) 
e. . 

E(X»t ,03) 

g e (y,a,t) 

g l (y,a,t) 

G s (Y,A,T) 

H ijU 

1 A’ 1 B 

*A’*B 

X o 

k 

K 

£ 

L 

M 

M 


P ij 


P ijki^- , - ,T ' ) 

£ 

r 


r e (Y,A,T) 


ambient sound speed 

acoustic intensity autocovariance, equation (3.6) 
rate of strain tensor 

acoustic energy density, equation (3.6.1) 

Eulerian corssed-beam covariance 
Lagrangian crossed-beam covariance 
Schlieren covariance 

spectral density transform of equation (3.18) 

detector signal fluctuations, beams A and B 

mean detector signals, beams A and B 

detector signal , no extinction 

fluctuating extinction coefficient 

mean extinction coefficient 

distance variable 

turbulent scale length 

convection Mach number vector 

IMl 

pressure 
stress tensor 

Lagrangian scalar covariance 
Lagrangian stress tensor covariance 
vector distance from disturbance to observer 
111 

component of r 

Eulerian scalar density covariance 
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R ijki^— ,T ^ 
S(Y,t,T*) 


t 


T ij 

u i 


u 


U 


U 


v 


V, 


V 


r 


V 


X 

y s 

Y 

Z 


A 




e 


K 


X 


r 


Eulerian stress tensor covariance 

sound source intensity 

time 

Lighthill stress tensor 
mean velocity 

velocity component in observer's direction 

mean convection speed 

convection speed 

velocity vector 

component of v 

component of v in direction r 
turbulence volume 
component of observer's position 
vector locating observer 
component of turbulence position 
component of Y in observer's direction 
component of Y normal to r 
vector locating turbulence source 
vector locating turbulence source 
frequency parameter 

(Y - Z) separation between points in turbulence 
component of A normal to flow direction 
component of A normal to 
angle between M and r 
wave number vector 

Lagrangian distance variable in r direction 
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A 


^n 


A 

— s 


Po 

a 


x 

T* 

T 

CD 


Lagrangian separation between points in turbulence 
component of A, normal to M 
component of A normal to r 
density 

ambient density 

Lagrangian time variable 

retarded time difference 

arbitrary time delay in sound source 

arbitrary time delay in crossed-beam signals 

frequency 


Operators 

V2 

□ 2 

f 

( ) 


Laplacian 
D f Alembertian 

indicates equality of two functions after integration 
time average 
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APPLICATION OF CROSS ED -BEAM TECHNOLOGY 
TO DIRECT MEASUREMENTS OF SOUND SOURCES IN TURBULENT JETS 


I. INTRODUCTION 

(l 0) 

The crossed-beam technique has been shown' * y to per- 
form an area integration of correlations very similar to inte- 
grals appearing in the source term of the equations for turbulent 
noise production. This suggests then, that there is a possi- 
bility of applying crossed-beam techniques to the direct measure- 
ment of local sound source intensities in turbulent jets. The 
advantages of being able to do so are obvious : it becomes 

possible to study the basic nature of noise production and the 
effects of various noise reduction schemes upon the source itself. 
Indirect inference of sound source characteristics from far-field 
acoustic measurements can never rest on very firm ground. 

This report represents a study of the mathematical con- 
cepts involved in applying crossed-beam technology to source 
measurements in order to indicate the conditions under which 
such an application is possible. For this purpose, much of this 
report is devoted to a redevelopment of basic concepts conceived 
by others. This is necessary to give the reader the required 
background and to present the pertinent equations in forms ac- 
cessible to crossed-beam comparisons. The basic crossed-beam 
equations have been developed in several places, and are re- 
viewed briefly in the following section. The approach to noise 
generation parallels quite closely to Ffowcs-Williams ^ four- 
dimensional Fourier transform approach since the equations in 
that form are manipulated readily for our purpose and point up 
the physical characteristics clearly. 
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II. CROSSED-BEAM PRINCIPLE 


The operation of the cross ed-beatn technique has been 
described elsewhere and will therefore be discussed only 
briefly here. In principle, two narrow beams of radiation at 
selected wavelengths are positioned such that they intersect at 
a point in the region of interest in the jet. If the beam in- 
tensity of each beam is modulated by a process which is linearly 
related to the local gas density fluctuations, then the inten- 
sities can be expressed, apart from a calibration constant, by 


i 


A 



( 2 . 1 ) 


'B 



( 2 . 2 ) 


where i^, ig, 1^ and Ig are the fluctuating and mean signals, 
respectively, for each of the two beams A and B. If the two 
signals are cross-correlated over a sufficiently long inte- 
gration time, the product mean value ^i^ ig^> will give infor- 
mation concerning only the fluctuations common to each beam 
(i.e., in the intersection region). In this way, local infor- 
mation on turbulence has been obtained without the necessity 
of inserting physical probes into the flow field. From a more 
intuitive viewpoint, we might consider turbulence as consisting 
of discrete eddies that are slightly more or less transparent 
to the beams of radiation than the surroundings. As an eddy 
passes through a beam of radiation, the light intensity at the 
detector changes. When the two detector outputs are cross- 
correlated, only the information concerning eddies passing 
through the beam intersection point is retained. Signal fluc- 
tuations due to source instabilities, detector noise, etc., 
are normally uncorrelated, although these contributions do 
affect the overall signal-to-noise ratio of the system. 


NT RESEARCH INSTITUTE 


2 



Cross-correlation of equations (2.1) and (2.2) gives 


G e (Y, T) = 


< c i a 

t a t b 



y 2 + *2» y 3’ t ' ) p ^ y l 


+ y 2 , y 3 + V* + T >^ dA 2 dA 3 (2.3) 

where i A and i g are correlated with an arbitrary time delay 
between signals, and A^ is a displacement of beam B relative 
to beam A in the direction of their common normal. This equa- 
tion expresses the correlation between two line integrals. If 
we can write equation (2.3) as 


G E (Y, A x , t) =JJ<^p(y 1 , y 2 > y 3 > t) P (y x y 2 +& 2J y 3 


+ A 3 , t + T)^> dA 2 dA. 


■/ 


r e (2. t) d^ 

where d^ = dA2 dAo or, equivalently as 


(2.4) 


G e (Y, a, T) = <(p(Y, t) j p(Y + A, t + T)^> dA 


— tl 


(2.5) 


we are then dealing with the correlation between the area inte- 
gral of points in a plane (y^ + A-^ = constant) with an outside 
point defined by (y^, y 2 » y 3 ) • Such a cross-correlation has been 
assigned the term "point-area correlation and is similar to 
the type occurring in turbulent sound source integrals . 
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If the turbulent flow is everywhere homogeneous , then 
the point-area formulation is obviously proper. If departures 
from homogeneity are pure odd functions of space coordinates 
about the point (y x , y 2 , y 3 ) , then these departures tend to 
integrate out and the point-area concept is still validP^ In 
any event, the final proof of the point-area correlation assump- 
tion can come only from experiment. These experiments are the 
subject of a separate study. 
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III. THE EQUATIONS GOVERNING TURBULENT NOISE PRODUCTION 

LighthilL^ has developed the theory of aerodynamic 
noise by the use of an acoustical analogy, rearranging the exact 
equations of motion into the form of the inhomogeneous wave 
equation with a source term. If there are no sources of mass 
or external forces present, the source term becomes quadrupole 
in nature, existing in the form of the double divergence of 
momentum flux. 


The wave equation becomes 


□ 2 


.2 = 


p a? ' a ° ' p 


° X ij 


where 


: ij = P v i v j + P ij ' a o P 5 


3-J 


(3.1) 


(3.2) 


from the exact equations of motion and represents the quadrupole 
strength. This equation can be written 


(p - p D )(X, t) 


i r d \i <*•* 

4-n a c J dy^y^ 




(3.3) 


Equation (3.3) is actually only the integral form of equation 
(3.1), not a true solution, since the stress tensor T. . is it- 
self a function of p. The only way a solution can be obtained 
is to approximate T^ by terms independent of the density ap- 
pearing in □ 2 P . The foremost difficulty in using equation (3.3) 
for noise calculations , is in fact choosing the approximation 
for T i j appropriate to the problem of interest. 

The "solution" equation (3.3) is given with the inte- 
grand to be evaluated at the retarded time |X - Y|/a which is 
the time difference between time of emission at Y and time the 
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sound wave reaches the observer at X. This retarded time ef- 
fect is the major stumbling block in casting the noise equation 
in a form identical with the crossed-beam equation. However, 
it must be pointed out that if retarded-time effects did not 
exist, equation (3.3) would integrate directly to zero for the 
far field (large |X - Yj) as will be shown later. Physically, 
this says that the quadrupole (or dipole) type of sound 
source creates a far-field disturbance only because of the 
phase mismatch which results from the varying distances to 
individual disturbances. This is the Stokes effect of classical 
acoustics . 

Equation (3.3) can be reformulated in a more convenient 
way, and one which brings out the above points more clearly by 
separating out the retarded time variable using the definition 
of a delta function, i.e., 


+CO 


(p - p o )(x, t) = 


<^T y Qf, t) 6 (t - Xj) 


4tt 


i f f S T ii<: 

—1 I dY dt 

a °-u„ y i y j 


where 


X - Y 


a. 


_ r_ 
a. 


is a function of Y. 


Rearranging terms , 


(p - p Q )(X, t) - 


*5 // ^ 

■If 


St. . 


r dyj 


^ dY dt 


St . . % t 

'SyT d ~ dt 

3 J J 1- 


(3.4) 
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Now the first term can be converted to a surface integral of 


6 

r 


St. . 

377 1 


by the divergence theorem. If there are no surfaces present in 
the flow, this must vanish at infinity. 

If we repeat this process we get 

(P - p 0 ><£. t) ‘ T ij ftsyj <£> ^ 

Differentiating by parts 

-s c 6r . . 6 1 r . 

a /6 k 1 1 1 

5y 7 ^r^ 3 a " "T ’ 

■'l r or 


where 6' = d6/dr. Taking S/Sy^ of this and substituting into 
the above equation, we obtain 


(p - p o )(X, t) 


-V f 

h'na Q J St 



i r 

7 


dY 


2 

where we have dropped terms of order 1/r and less (near field 
terms). This may be expressed as 


(P " P G ) CX , t) 


7 77 

o 



(Xj - y t ) (x 1 - y.) 

[X - Y | 3 



(3.5) 

The autocovariance of the far field density fluctuations is then 
obtained by '’squaring" equation (3.5) and averaging in time: 
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<^P - P c > (X, t)(p - p 0 )(X, t + T*)^ - B(X, t, t*) 

= i f f < x j-yp< x i ~ yi > < x k ~ z k> -*i> v 

16**8® J J _ |x - Y| 3 |X - Z| 3 

✓ *ki<X.V * k i<&.* 2 + **>V ... 

/ dl « 


(3.6) 


x* is an arbitrary time delay, introduced only to be able to 
express B as an autocovariance. The acoustic energy density is 
then obtained from the Fourier transform of B 


E(X , t, o>) 
Substituting 

T = ^2 + T* 



t , 


T *) e -i“ T * 


dx* 


X - Y| - IX - Z 

X, = T* + 

1 a_ 


(3.6.1) 


A * (X - Y) 

= t* + (3.6.2) 

- H 


and assuming A ^ |X - Yj in equation (3.6) 
B(X » t, x*) = 


i fr (x i- y i )(x i v 

T^JJ ^ x 


A, T) 


dA dY 


(3.7) 


where 


R ljW a. A, t) = <T^j (Y > T X ) T W (Y + A, T ± + r)> (3.7.1) 

and is to be evaluated at the retarded time x. x is a function 
of A and is, consequently a dependent variable. It must be 
treated as such in any manipulation of equation (3.7). 
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In a turbulent flow, the radiating quadrupoles are not at 
rest and move during the retarded time t . This motion is cor- 
rectly accounted for in equation (3.7) since the correlation 
tensor measured in the frame of the observer includes the 

convection effects. If we consider the time derivative of such 
a correlation, we find there are two contributors, the true time 
change and the change due to convection of space derivatives. 

In most turbulent flows, it is the latter which dominates, i.e., 
apparent changes with time are due to the instantaneous convec- 
tion of space derivatives. Even though these space derivatives 
may be very large they do not contribute to the integral of equa- 
tion (3.7). This can be proven in the fashion discussed earlier 
by expressing volume integrals of divergence terms as area inte- 
grals which must vanish for large A if 

lim R iikX^— ’ — ’ ^ ~ 0. 

Thus, in order to present equation (3.7) In a form use- 
ful for noise prediction we must write it such that pure time 
effects are maximized. This can be done by introducing a moving 
frame of reference in which time scales are maximized. This of 
course is the frame of reference moving at the convection speed 
of the covariance. It must be stressed that this only minimizes 
the effects of convection of derivatives. They can be zero only 
in a frame which moves instantaneously with the fluid. The con- 
tribution of the velocity fluctuations poses another problem to 
us in the measurement of sound source strengths. 

We note here that if a reference frame exists in which 
= 0 

then there can be no sound. This means that a convected field 
of "frozen" turbulence does hot radiate sound. There are some 
very important exceptions to this and these will be discussed 
later. 
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Transforming to a moving axis using 

* = A + a G M t (3.8) 


where M is the convection Mach number measured relative to X, 
we define a moving-axis correlation tensor by 

P ijki<X> ^ t) = R ijki ( £’ T) (3.9) 

The retarded time (see sketch of Figure 1) is then given by 

(X - Y) 

a 0 T - a o x* + (A + a M x) • (3.10) 

is - 1| 

or 

A • (x - Y) + a o T* JX - Y| 

a o {(S - I| - M • (x - Y)^ 

The time derivatives are related by 


^ R iiki! _ fd d 1 

3t jot ' a o M r> 35^ | p ijW (3.11) 

Recall here, that t is a dependent variable (t = x(?\)) 
and so d/d>^ includes some time effects; it is not entirely a space 
gradient. These effects are determined using equation (3.10) 
as follows: 



p ijki<X, A, t(A» 




P ijki! (I> I’ T > 


d_ + Sa " y n 

a 0 [|I ’ I| - M ‘ (X - Y)] 



t) 


(3.12) 
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X 



Figure 1 COORDINATE AXES SYSTEM 
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Here, P. . , (Y, X, t(X)) signifies that t is to be considered 

a function of X prior to differentiation and (X> X» t ) 

signifies that t is to be considered an independent variable 
until after differentiation. Using equations (3-11) and (3-12) 
we get the final relation between time derivatives in the 
fixed and moving frames 


_S_ D 

St ijk£ 


(I, A, x) = 


|X - Y 


X - Y | - M- (X - Y) 


_S_ 

Sx 


- a M 


o n SX, 


h 




(3.13) 


Again, the space derivative 


SX, 


P. 


(X) will contribute 


- n ijU ^ 

nothing to the integral of equation (3.7) and we will drop that 
term. Differentiating three more times, equation (3.13) becomes 

[4 A 


_S 

S T 


R ijki 


(X> A, t) f 


IX = Y 


[IX - Y | - M- (X - Y)j St 


P ijkJ5 


(X, t) 

(3.14) 


where f signifies that the two sides of the equation are 
equivalent when integrated as in equation (3.7). This is, both 
sides of equation (3.14) produce the same far field effect 
although they are not strictly equal. 


The Jacobian of the transformation gives 

IX - XI 

dA = dX (3.15) 

lx - xi - x- (x - X) 
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Substituting equations (3.14) and 3.15) into (3.7) gives 
finally 


B(X, t, T*> = 


16" 


1 f <* -s'i>(*i-yiH*k-yk> (x e- y e ) y 
^ J |x - V|[|x - Y|- M-(X - Y)] 5 



P ljla <Y, i. T) dX dY 


(3.16) 


It is the inner integral which will concern us here since it 
contains all the local source characteristics. We will call 
this integral 


S(Y, t, t*) 



P ljW (Y, i, T) dX 


(3.17) 


The characteristics of this source integral have been dis- 
cussed by Ffowcs -Williams by Fourier transforming S(Y, t, t*) 
to obtain its sound power spectral density tensor’. In this 
way we can see more clearly the effects of source frequencies 
and physical scales on the resulting sound frequencies. 


The four -dimensional Fourier transform of the correlation 
tensor, defining a power spectral density, is given by 


H ijki (I, K, ») = 

from which 
P 


M 


p ijki ° )e 


•io:a , 


X e * — dX do 


(3.18) 


*ijki ’ -» o) ^ J H ijk£ 


(Y , k_, co)e e — — die da> 


(3.19) 
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The four-dimensional transforms are used since we must 
determine the effects of both wave number (space) and frequency 
(time) effects within the source and the relation between them 
under the restraint of the retarded time relation, equation 
(3-13). In the far field a simple dispersion relation exists 
between k and cd and a one-dimensional transform will suffice. 


Differentiation of equation (3-19) with respect to o is 
equivalent to multiplication by ico so that we can write, using 
equations (3.10), (3.17), 3.19) 


SQL, t, t*) = ( 2") 3 I I go 4 H. jki , (Y, k, cd)X 


/**, 

( 2") 3 J j 


(COT + £ • X) dx dK 


CO 


H ijkjg QL> a >) exp 


ico IX - Ylx- 


X - Y| - M* (X - Y) 



(2”) 3 



co(X - Y) • X . . 1 

a o [l— - II- M-(X-X)]j 


dK 


dco 


H ijk^ 


(I, 


-co (X - Y) 

a Q [l— * II - M- (X - Y)] 


co) 


x exp 


r ico|X - Y|t* 

; |X - Y|+ M- (X - Y) 


dco 


( 3 . 20 ) 


As an interesting aside, the calculation procedure here is 
similar to that we do to reduce the crossed-beam power 
spectral density to a one-dimensional Fourier transform, 
the difference being that the crossed beam involves an area 
integral of a Fourier transform, whereas here we have a 
volume integral. 
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At very low speeds, M — ► 0 at a Q -^a> and = 0. 

Thus at low speeds for fairly low frequencies ^~0 ° 

by equation (3.20.1). The frame of reference chosen ensures 
that time scales are maximized and hence frequencies will be 
limited. 


If we take k = 0, then equation (3,20) gives 
3 


S(Y, t, t*) - (2") 


) f 4 fH£ - il T * 1 

I “ (x> °* m) Hx-Il-M-(X-Y)/ 


do) 


(3.21) 




P ijk£ 


T*|X - Y 

Y, A, 

X - Y | - M* (X - Y) 


dX 


(3.22) 


From these equations we can see several important charac- 
teristics of turbulent sources at low speeds. 

The frequency is increased by the factor 

|X - Y | 1 

|X - Y | - Mi (X - Y) 1 - Mcos 0 

a result of motion of the sources relative to the observer. 

This is just a doppler effect. 

The acoustic energy at low speeds is independent of the 
retarded time t, which was precisely Lighthill's low speed 
assumption. The sound sources are compact in nature in 

the moving frame of reference so that retarded time differences 
are small; we have accounted for retarded time differences 
resulting from convection by the coordinate transformation. 

This is equivalent to assuming that all the acoustic power 
comes from zero wave number turbulence or, what is equivalent, 
the power comes only from true time varying signals (in the 

moving frame of reference) . We must note here that k ~ 0 is 
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an approximation only. In actual fact if k = 0, there would 
be no sound since co = 0 then too. 


As the convection speed increases, k can no longer be 
neglected since [X - Y J - M. (X - Y) decreases and in face a 
singularity occurs in equations (3.20) and (3.20.1) when 
|X - Y| - M* (X - Y ) =0. The occurrence of the singularity 
introduces an interesting and important aspect of turbulent 
noise. Even though equation (3.16) appears to indicate that 
frozen turbulence creates no sound, we can now see an excep- 
tion mentioned earlier. The singularity results in an 
indeterminate of the form 0/0 which can and does have a 
finite limiting value. The singularity condition is the 
case where the turbulence convection speed has a component 
toward the observer equal to the speed of sound. This is 
precisely the condition required for sustaining a Mach wave. 
Thus the eddies are analogous to solid bodies creating Mach 
waves which propagate toward the observer, and hence the 
name Mach wave emission for this type of (supersonic) turbu- 
lence-generated sound. 


The retarded time difference becomes infinite at 
the singularity and each source component of the quadrupole 
then behaves as a separate source, no '’Strokes effect" can- 
cellations occurring. Thus Mach wave emission is very effi- 
cient (sourcelike). It will radiate for all M = — ) 1 at 
the Mach angle 9 relative to the flow direction ° 
defined by 


M- (X - Y) 


X 


Y 


= Mcos 0 


(3.23) 


The simplified expression (3.22) obviously does not hold 
for high speeds and we must go back to equation (3.20) which 
holds for all M. 
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Introducing a new frequency variable a into equation 
(3.20), defined by 


cu ] X - Y | 

a = 

X - Y | - M* (X - Y) 


(3.24) 


we get 


S(Y, t, t*) 



X - Y| - M- (X - Y) 

|X - Y| 


5 


X H ijk 
X e -iaT * 


a(X - Y) 

V 

X - Y -M- (X - Y) \ 

— 5 1 , 9 

a oi*-H 

|X - Y ) 


da (3.25) 


The singularity in equation (3.6) is cancelled by the corre 
sponding term here and near the singularity the frequency term 
becomes zero so that 


H ijk^ H ijkJ 


I, 


-a (X - Y) 

,0 

a ol^- H 


(3.26) 


This indicates that it is the wave number K r in the 
direction of propagation 


(X - Y) 

il - II 

to the observer along the Mach acoustic ray at zero frequency, 
which produces the sound power. Note that cd = 0 is required in 
order that a from equation (3.24) have a finite value. It is 
the physical scale rather than the time scale of the eddy that 
becomes important, the complete antithesis of the low speed 
case . 
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Fourier transforming equation (3.25), 



where X r is the component of X in the direction of the observer 
corresponding to the wave number component k = — in that 

IT £L 

direction. ° 


Then by equation (3.25) 


S (Y , t, t*) = (1 - Mcos 0) J -°Mo 


5 a) 


X 6 


T* + 


X- (X - Y) 


a cl* - H , 


da dX 



<y P • 

= (1 - Mcos 0) 5 a^ j j—2-LiM (y, x, a) do dX g 




where X is in the area normal to X defined by 


(3.28) 


(X - Y). X s = 0 
and where 


(3.29) 


X r * - a o T * 


(3.29.1) 


The important frame of reference then is one with two 
axes in the plane tangent to the Mach front and moving with 
it . 
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Thus as an eddy moves it emits continually as long 
as there exist a space-derivative fluctuation. It will 
continue to emit throughout its life-time which is finite, 
otherwise the time integral would give infinite power. The 
lifetime here is the timescale of The smaller are the 

changes with time obviously the larger is the acoustic power 
radiated, in marked contrast with the low speed quadrupole 
case where changes with time are essential to the radiation. 

Even at supersonic speeds there also exists quadrupole- 
type radiation at angles other than cos'^ . Since however, 
the quadrupole efficiency is so much smaller than that of the 
source- like Mach emission, the former may be ignored for the 
most part, except perhaps at far downstream observer positions. 

In review, we have been able to write two integral 
expressions where retarded time differences can be accounted 
for readily. This ability is essential to us in applying a 
space- integrating device like the crossed beam to sound source 
measurements since different parts along each beam cannot 
contribute at different retarded times. The two expressions 
are applicable to 

1. low subsonic speeds, where retarded times can be neglected 
in a frame of reference moving with the eddy and 

2. high speed flows in which the ambient sound speed is 
somewhere exceeded by the eddies. 

In this latter case the proper frame of reference is 
one in the plane of the Mach wave front and moving with it 
at the ambient sound speed toward the observer. Since points 
on this plane are equidistant from a far field observer, 
retarded time differences vary negligibly in the plane. Since 
the eddy (and the reference frame) travels at the speed of 
sound toward the observer, Stokes cancellation effects and 
retarded time changes in that direction do not enter the 
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formulation. In both the low speed and the high speed case, 
the correlation tensors must be measured in the moving ref- 
erence frame specified. 


The equations for source intensity are: 
High Speed -- Mach Wave Emission 


S (Y , t, t*) = (1 - Mcos 9) 5 a^ 


Low Speed - M-^0 


■Hg* 


(X» A*. o) 


dX d o 
— s 


( 3 . 28 ) 


S(Y, t, t*) 


P ij ki 


X, A, 


t*|X - Y 


dX 


IX ‘ xi - X* (X " X) 

( 3 . 22 ) 

Under other conditions we must use the general expression 

,4 


S(Y, t, t*> 


-/ 


cK 


- (21T) I “ H ijkg(I> 


P ijM (X, A> T > d A 


- 03 (X - Y) 


( 3 . 17 ) 


a o[|X * X| - M. (X - Y)] 


> /J 0 X 


exp 


ioD |X - Y | t* 1 

. X " XI - M- (X - X> J 


doj 


where 

X* (X - Y) + a Q t* |x - Y | 
a 

o 

-03 

” a o[lX - XI - M-(X - Y)] 


[lx - XI - X* (X - X) ] 
(x - X) 


( 3 . 20 ) 


( 3 . 10 ) 


( 3 . 20 . 1 ) 
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and 


p ijki ( i’ i> T) 
(3.14) 

The arguments given for cancellation of retarded time 
effects with reference axes in frames normal to the observer 
would also be true in the general case of equation (3.17). 
However, a volume integral is of importance here and retarded 
times are important in the observers direction. Whether or 
not these can ever be accounted for are left to the following 
sections for discussion, along with general discussions of 
the applications of all these equations to measurement by 
crossed-beam techniques. 
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IV. APPLICATION OF CROSSED BEAM TECHNOLOGY 


The problems involved in applying crossed-beam technol- 
ogy to sound source measurements lie in three areas. 

1. Covariance Measurements — approximation of 
the covariance by parameters measurable by 
crossed-beam techniques. 

2. Retarded Time Effects. 

3. Moving Frame of Reference. 

These are discussed in this section with subcategories 
for delineation of (1) low speed applications, (2) high speed 
applications and (3) general applications. 


4.1 COVARIANCE MEASUREMENTS 


The signal measured by a single beam of the crossed- 
beam system is determined by the extinction coefficient of the 
dominant light extinction process. The detector signal I is 
given by 



If the intensity is assumed to result from a mean I and 
a fluctuation i, 


i + I = I Q exp 


A 


= I exp { 
i = T ( exp { 


(K+k) di 
kdi 


•kdi | - 1 

and if the total integrated fluctuations are small 
i(t) = I j I k(Y,t)di 


(4.2) 


With proper calibration constant accounted for, the single or 
crossed-beam system then measures the statistics of the extinc- 
tion coefficient whereas we wish to measure the statistics of 
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a fluid property. Within a gas, and for most extinction proc- 
esses available to us, the extinction coefficient is strongly 
dependent upon the number of contributing particles or to the 
mass density if the particles causing the extinction are uni- 
formly distributed on the average. The present generation of 
crossed-beam systems uses absorption by naturally occurring 
CC >2 in the infrared region of the e.m. spectrum for the extinc- 
tion process. Separate studies are underway to determine under 
what conditions of wave length, etc., k is proportional to p, 
and independent of pressure and temperature. 

In any event it appears that density statistics are the 
most promising approach for crossed beams and we must therefore 
be able to express the sound source covariances in terms of 
density . 

The source term^^ involves the quadrupole strength 

T ij ■ ? v i v j + P ij - a o 2 P«lj 

2 (4.3) 

~ P V i V j + (P “ a Q P)*ij 
if we neglect viscosity. 

In isothermal flows 

2 

p ~ a Q p, and Tjj ~ P 0 v i v j (4.4) 

If there are large entropy gradients in the flow, however, the 
pressure and density fluctuation may not cancel out and p-a p 
will be determined by entropy fluctuations. 

Thus in a heated jet equation (4.3) might become 

~ (P- a o 2 P> a ij < 4 - 5 > 

In the case of a supersonic jet, of course, all terms must be 
included in T^j . 
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4.1.1 Low Speed Jets (M=0) 

We will restrict our discussions to isothermal jets at 
low speeds. This is obviously not because of a lack of impor- 
tance of heated-jet noise but is rather only to limit ourselves 
to the simpler case. A much more comprehensive study is required 
for the entropy fluctuation case and this is left to a future 
time . 


Shear Flows 


In most turbulence flows of interest to us here, there 
will be a large shear gradient in the flow, and turbulence will 
tend to be most intense in that region. If shear is present it 
has an amplifying effect upon the noise produced. From the 
momentum equation 


at^Vj) = 


p ij 


p- 



(4.6) 


The last term will drop out when integrated over all space, 
(using the divergence theorem) and dropping viscous terms we 
get from equation (4.4) 



(4.7) 


The mean rate of strain tensor may be approximated by its mean 
value which we will assume to have a dominant direction y-^ 


e 


ij 



(4.8) 


Then 



(4.9) 
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(4.10) 



(4.11) 


Taking the mean-square value of this, 


a 


4 



R ijki 




) 


where 


R e (Y,A,t) = <p(Y,t) p(Y + A, t+r) > (4.11.1) 

We are therefore concerned with the second derivative 
of the density covariance. It must be expressed in the moving 
reference frame for use in equation (3.28), but this will be 
left for later. 

Homogeneous Turbulence 

If there is no appreciable shear in a flow we can ob- 
tain R. ... in terms of the density covariance using the approach 

1J K JO / q \ 

of Ribner.^' He considers the pressure field as composed of 
an ambient p , an incompressible part p and a compressible part 

o 2 

p with corresponding densities related by p = a Q p . 

For incompressible flow Dp/Dt = 0 and equation (3.1) be- 
9 9 

comes -V p = (d T ij /^ x i^ x j ) inc where T ij is evaluated at zero 
convection speed. Substituting into equation (3.1) we obtain 



V 2 p - 



V 2 p° 


(4.12) 


where we have assumed (T i j ) comp " ( T ij>i n comp at zer0 convec " 
tion speed and 


o 

P 



o 
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The source strength then becomes 



iikl 


f 


with 




R(Y,A,t) 


R(Y,A,t) = < p(Y s t)p(Y + A, t + t)> 


(4.13) 


(4.11.1) 


Thus at low speeds the source strength can be written in 
terms of the density covariance as required for crossed-beam 
measurements . 


4.1.2 High Speed Jets (M=l) 

From equation (3.28) the primary term in the sound source 
integral for Mach emission is 


with x r 

OA, 


Now by equations (3.4) and (3.16) 


constant 


(x i - y i )(x i - y.|)(x k - y k )(* { - y t ) 

r b (l - Mcos S)-" 


ax 4 


a 2 T.. (Y, 0 ,t) 

< ‘ L 

hy^dy^ 


a 2 T t1 (Y,X, 


t + t) 

> 


and we are concerned with the component of this in the r 
direction. 


Now, 


l 

r 


iin. 


f v -i 


r a 


t = 


r_ 

a 


2 

a T 


at^y i 


11 = 


j-ih 

ra Q btby t 


(4.14.1) 
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5T ir . 

W 


dP 


Sy t 


<P v i v j + P ri> ‘ a o — 

3p 


2 °P v r 

-a n H 

° ^y r St 


(4.14.2) 


We may write 


1 

r 


-1 


dpv r 


dy r 


a r St 
o 


(4.14.3) 


Substituting equations (4.14.2) and (4.14.3) into (4.14.1) we ob- 
tain 


S 2 T . . 
LI 

Sy^y. 


f ^ 


a_ 

St 


a & - ^ 

° 3y r Sy r 


r 0 f , Sp ^ V r 

* ■ * <a ° ■ v t ■ p * 


If we let v = u + v in equation (4.14.4) we get 


( 10 ) 


s 2 t. . 


Sy t Syj 


= (a - u ) 
v o r / 


n .2 

a P 


a 2 pVl 


Sp Su 


+ — 


Stdy r StSy r dt ay 


(4.14.4) 


(4.15) 


Shear Flows 

If we are dealing with a high-shear region, only the 

third term is amplified by that shear. The other terms will 

tend to be small. Particularly, we can see that with U -a 

co 

that the factor a Q - u will be small. Thus we will take 


d2l n _ au r ^ _ <*i -y 1 ) (x i -Yj) 9u t 9p 

Sy^y^ ay r at r 2 a yj at 


(4.16) 
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where we have substituted 


a u 


9y r 

This gives 


(*i - y i )(x 1 - y.) S Ui 

2 U 7 

r 


^ p i ik ^ (X»*. T > = 


a* H 


Su i\ / ^ p OC>t) Sp (X. + t + T) 

l ^ 


dt 


y (4.17) 


Since this is to be determined at A = const, and further, 

du. ^ du x r 

taking = equation (4.17) becomes 


dyj dy 2 


1 du x ' 2 


ax 4 s 


\ d ^2 


i?" 


A r = a o< T “ T *> 


(4.17.1) 


Thus 

S (Y 5 t*) = - a ^(l - Mcos 9) 5 


ft 


du ,\ 2 a 2 R F (Y,A t) 


dy 2l 


dx z 


dA g do (4.18) 


where 


R E<I>V T) = <p(I. t )p(l + 4, t +t)> 

A = a (t -t*) 
r o v x 


(4.18.1)- 


Except for the sign, equation (4.18) is by coincidence identi- 
cal with equation (4.11) for low-speed shear flow. The important 
observation here is that the source intensity is determined by 
the area integral of an Eulerian covariance. The gradient terms 
in the time derivative do not integrate out in an area integral 
and in fact dominate. 
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Homogeneous Turbulence 


If the flow of a jet is supersonic there must be 
large shear gradients in the flow as the velocity goes from 
zero outside the jet to greater than speed of sound within 
it. There should be no necessity then of attempting a 
homogeneous solution for the Mach emission case. 

4.1.3 General Case 

In the general case, when we cannot assume either 
M — 0 or M cos9 = 1 , it is not possible to approximate 
by a simple function of density. As stated earlier however, 
the Mach emission case probably covers all speeds of interest 
above M=l. Below M=l, the low-speed isothermal approximation 
is no doubt valid for isothermal jets over a broad range of 
speeds. More serious of course is the need for an expression 
in the case of heated flows. It appears that it should be 
possible, however, to express entropy fluctuations directly 
in terms of temperature, and consequently density fluctuations. 

In an isothermal jet, it is not so much the approxima- 
tion of the stress tensor covariance that presents the most 
severe problem. Rather, the retarded time differences 
become nonnegligible as convection speed increases , even in 
the moving frame of reference. Means of accounting for 
retarded time differences are discussed in the following 
paragraphs. 

4.2 RETARDED TIME EFFECTS 

In section II, it was shown that the crossed-beam 
system measures a covariance given by the expression 

g e (y,a 1 ,t) m JJ< ( p(yi •y 2 + A 2 ,y 3 + A l’ y 2’ y 3 + A 3 > c + T )> 

X dA 2 dA 3 
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Under conditions of homogeneity or odd-function departures 
from it , this can be written as a point-plane correlation 
leading to 


g e (y,a 1 ,t) 


r E -(X.£,T) dA n 


(2.4) 


where dA R = dA 2 dA^ is a plane normal to A^. The sound source 
integral on the other hand is generally of the form (from 
section III) 


-/ 


S(Y ,t ,r*) =| -^T P ijk dA 


£ (1-Mcos0)~* 


J dr 4 


t) 


dA 


(4.19) 


Comparison of equations (2.4) and (4.19) suggests that we might 
obtain S(Y,A,t*) from a crossed-beam measurement using a 
relation of the type 


S(Y,t,x*) 


(1 - M cos0) 




(X^^T) 



(4.20) 


We can see now the need for obtaining the crossed- 
beam covariance in the form of the point-plane covariance is 
to approximate the correct covariance R under the area 
integral. 

A difficulty involved in the use of equation (4.20) is 
immediately apparent; x must be made identical with T, so 
that 


T » T 


T * + A_cos0 


(3.7.1) 
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In order that Gg(Y ,A^ ,T) be expressed in its point- 
plane form of equation (2.4), it was necessary to assume that 
we could make arbitrary coordinate translations normal to y > 
of the type y 2 “ ^y£ “ ^2 w ^ t * lout changing the statistical prop- 
erties of the covariances . This is a good approximation prob- 
ably, but when the translation is made we also vary x* through 
equation (3.7.1). It therefore becomes impossible in general 
to maintain the correct x for arbitrary 9 and ^ . There are, 
however some special cases where we do so and, these are dis- 
cussed in the following paragraphs. 

Another problem is associated with taking the time delay 
derivative of Gg. In equation (4.20) we have assumed that 


/ 


a n R E (Y,A,T) 


3T 


.n 


dA 


a n r 

SF J R ^ 


T) dA 

— n 


,n 


st 


- g e (y,a 1 ,t) 


(4.21) 


This is obviously not the case as long as T = T(A) which in 
general is the case as required by equation (3.7.1). Again 
there are special cases where this is done and these are 
discussed later. 


4.2.1 Low Speed Jets 

As shown in section III, if we transform to the moving 
frame of reference a large portion of the x differences are 
accounted for by the convection equation. The time scales 
in the moving frame become very large in comparison with the 
retarded time differences ~L/a Q , where L is a typical correla- 
tion length. If a G »U, then retarded time differences may be 
neglected in comparison with fluctuation time scales. Flows 
which satisfy this condition have been said to contain "compact" 
sound sources In this case then (U ^ a Q ) we can put x = 0 

and equation (4.20) becomes valid. 
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Shear Flows 


For high-shear flows , 


d 4 R 

17 


ijkl 


f 



7r e (y,a,t) 

77~ 


(4.11) 


or in the moving reference frame 


or 



ijkl 


f 


du. , d 2 P(Y,A,T) 

dy 2 a ° dr 2 


S(Y,t,r*) 


4 / du i \ 2 f ^( 1 , 4 , t) 

a ° tel J 7 ? 


dA 


(4.22) 


If G(Y,A-^,T) is measured in a frame of reference moving 
with the convection speed of the flow and the flow direction 
is along y^, then integral of equation (2.4) can be evaluated 
T = t = 0 . Then T ^ T(A) and we can use equation (4.24) 
to give 


S(Y,t,x*) 



(X.A-pT) dA x 


(4.23) 


The means of determining this moving axis covariance 
is treated in a later section. 

Homogeneous Turbulence 

For homogeneous turbulence we use equations (4.13) and 
(4.21) to give 

S(Y,t,T*) = G l (Y,A 1 ,T) dA 1 (4.24) 

Again G^ must be measured in a moving axis frame. 
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4.2.2 High Speed Jets 

In section III the source integral for Mach wave 
emission was developed to be 


S (Y , t ,t*) - (1 - M cos0) D a; 


5 ff &4p i 

°JJ ^ 




d *s 


d a 


(3.28) 


where A is tangent to the plane of the Mach front and A 
is the distance variable in the direction of the ray. In 
this form the retarded time is approximately constant for 
the plane A r = const, and only enters into the problem for 
the time integration. Thus the surface integral itself has 
instantaneously no retarded time effects making it of suitable 
form for measurement by crossed-beam techniques. 

In the case of shear flow, and the only one we shall 
consider here, equation (3.28) can be written using equation 
(4.18) as 


S (Y ,t ,x*) - -(1 - M cos9) 




2 S 2 R e (I.^ s .t) 

ST" 


dA da 

—s 


Assuming du^/dy 2 varies slowly with A r we can write 


S(Y,t,r*) = -(1 - M cos0) 5 a^ 


From equation (3 . 8) 


2 rr £> 2 r e (y,a s ,t) 


a// 


d t 


dA da 

— s 


(4.25) 


A = A + a t 
r r o 


(4.26) 


and from equation (3.29.1) 


A 


r 


-a o T* 
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Consequently 


(4.27) 


The retarded time difference then is constant along 
planes of A r = constant. 

Now if we align the crossed beams such that they lie 
perpendicular to r and are spaced a distance A r apart 


G E (X,A r ,T) = 


R e (Y,A s ,T) dA £ 


(4.28) 


Apart from the differentiation required, it is necessary 
to measure G such that T = A r /a o + t*. This satisfies both 
the requirement for time retardation as well as following 
the fluid motion at the convection speed toward the observer, 
i.e., the measurement is then made in the moving frame of 
reference. 

The time variable o in equation (4.24) is the true time 
in the moving axis system. It also satisfies the retarded 
time relation. However since retarded time is independent of 
A, s j the differentiation may come outside of the area integral 
in equations (4.24), (4.27) and (4.28). Thus we can write 


S(Y,t,r*) 


-(1 - M cosG) 5 a 5 

o 


S G E (Y,A r ,T) 

St 2 


T-t* = — 


(4.29) 


4.2.3 General Case 

The success in choosing the coordinate system (A ,A ) 
to minimize retarded time effects in the supersonic , Mach 
emission case, leads us to consider the possibility of achiev- 
ing like advantages in the general case. 
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To this end consider the coordinate system (y ,y ) 

IT S 

with y r in the direction of the observer. Then r^ = r^ = r 
and equation (3.16) becomes 


B(X,t, T *) 



5 S (Y ,t ,t*) dY 

(1 - M cos9)* > 


(4.30) 


with 



S(Y,t,T*) 



(4.31) 


where the inner integral now becomes similar to the required 
cross-beam area integral. We note here that t can be taken 
as constant within planes of = constant and therefore write 
T = T (\) or t = t(A^). By equations (3.6.2) 


T 


T* + 



Since P r is to be measured in the moving frame of 
reference in order to minimize density gradient effects, we must 
restrict the time delay of the crossed beam to 

A 

1 = 

a„ M cos9 
o 
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Letting x* - 0 , since it has no relevance for our analysis, 


T _ t = -£■ (M cosQ - 1) 

a Q M cos© ~~ 


(4.32) 


This is zero, only when M cos9 - 1 or a Q = <» . These are 
obviously the two limiting cases already discussed of Mach 
emission and very low speeds. 

We are left then with two alternatives at this point, 
either (a) let the delay be T 5 ^r^ a o and c ^ ec k to make cer- 
tain that gradient fluctuations do not dominate or (b) let 
T = A r^ a o M cos0 and assurae total error due to wrong retarded 
time is small. Alternative (a) is probably unacceptable if 
M cosQ is much different from 1 but might suffice for M cos0 
close to 1, whereas alternative (b) may extend the range of 
applicability of the low speed case to higher speeds. 

Retarded time errors which arise under alternative 
(b) above also prohibit taking the time delay derivatives out 
of the integral S. Thus we are not justified in taking the 
curvature of G^. It would be necessary to differentiate 
signals before correlating the two beam signals. 

In essence we are just concluding that some accuracy 
improvement in measuring sound sources in subsonic jets can 
be obtained by orienting the beams such that they are perpen- 
dicular to the direction of the observer and using equation 
(4.31). 

4.3 MOVING FRAME OF REFERENCE 

The covariances pertinent to the determination of 
sound source intensities must be measured in, or at least 
related to, a moving frame of reference given by the trans- 
formation 


A — A a Q M t 
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This is required in order to minimize convection effects on 
the measured covariances . We are faced with two problems 
in this regard. First, the inference of a Lagrangian co- 
variance from a measured Eulerian covariance and secondly, 
assuring ourselves that the minimized gradients are suffi- 
ciently small to be negligible relative to the "true" time 
changes . 

4.3.1 Minimizing Gradients Contribution 

are concerned primarily with two components of 

Tv 

Tvi = - a Q M t (4.33) 

A^ is in the flow direction. 

Tv r = A r - a Q m[cos0]t (4.34) 

A r is in the direction of the observer. 

The first of these corresponds to the reference axis 
most suitable for low speed jets and the second for the Mach 
emission and "general" cases. 

The covariances measured in the convection reference 
frame most certainly maximize the average time scales. How- 
ever, it is believed that a measured convection speed is 
determined primarily by the energy-bearing eddies, other eddies 
being convected at different speeds. The time scale is there- 
fore not necessarily maximized for all eddies , but only on 
the average. Those eddies which produce the sound (either k = 0) , 
or o> = 0 in the two limiting cases) are not necessarily the 

energy -containing ones. If they are convected at a speed U , 

s 

then 

[it] true ' It + U s ly[ (4 - 35) 

NT RESEARCH INSTITUTE 


We 

the vector 

1 . 

2 . 


37 



instead of 


[ 3t] measured “ It + U c 


(4.36) 


as is measured at the mean convection speed. Hence, we have 
a contribution of (U c - U g ) dp/d y;L which may exceed the true 
time variations and our measurement may be meaningless to 
sound source estimates. It is necessary therefore to deter- 
mine the contribution of these fluctuating gradients experi- 
mentally. We can show that equation (4.36) can be used to 
give 

_ ^ 2r E ( £’ A 1’ t) 2 ^Re ( - ,A 1 ,T) 

L ■ 9 ^ u T 




(4.37) 


- 2 2 

the error in o P/3 t 

is 

o (2) P - (U 2 ■ 

= u 2 ) 

r^\(Y,A 1>T ) 

dA^ 

j dA, 2 

— n 


d 2 R. 


Sa. 


(4.38) 


can be measured using crossed-Schlieren measurements^ 2 -^ and 
the result then compared with the measurement of 

da 2 • 

With proper data analysis the same assessment may be made for 

a \ 

i7" • 

The curvature of the Lagrangian autocovariance d 2 P/dT 2 
has been shown to be somewhat sensitive to signal filtering. 
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In particular, hot-wire measurements indicate a different 
lifetime at high frequencies than at low frequencies. There 
appeared to be very little sensitivity to frequency at low 
frequencies. However, crossed-beam spectra tend to peak at 
high frequencies in contrast to hot-wire spectra which show 
no peak, having a large low frequency content. It may be 
then that crossed-beam Lagrangian autocovariances are sensi- 
tive to band pass of the associated electronics. This must 
be assessed experimentally. 


4.3.2 Measurement of Lagrangian Covariance Integrals 
4. 3. 2.1 Low Speed Jets 
Shear Flows 

By equation (4.23) 


S(Y,t,0)= a* 



dA 


1 


(4.39) 


With the crossed-beam system, we normally measure 
G(Y,A^ = const, T) with as a parameter as illustrated in 
Figure 2. 

The envelope of this curve is G^Y,?^ = 0 ,T) = 

Gg(Y,A^ = U c T,T) the Lagrangian autocovariance. The convec- 
tion speed is then determined by A-^ = A^ = U c T = 0 evaluated 
at points on the envelope. 


In low speed shear flews we are interested in G^(Y,A^, 
T = 0) or rather in its second derivative with respect to 
time delay. At T = 0, A-^ = A^ and 


S(Y,t,0) = a£ 


/ du j\ 2 f d 2 G(YA 1} T 

te) i 



(4.40) 
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For purposes of evaluating this integral we can let dA-^ = U c dT 
by Taylor's Hypothesis and equation (4.40) can be written 


S(Y,t,0) = a* 




^<^( 1 , 0 , 0 ) 



dT 


(4.41) 


We have used the time variable a to indicate Lagrangian 
time where T refers to Eulerian. The integrand 


d 2 G L (Y,0 L ,0) 



is obtained from the curvature of G l (Y,}^, 0) in the sketch. In 
general G^(Y,A^,T) varies slowly and its curvature is very small 
necessitating very accurate measurements at each value of . 
There are always experimental variations in the heights of the 
covariances at each value of A^ and these may make derivative 
measurements impossible. Another approach is to differentiate 
the signals before correlating, to give 

^ g f(I> a i > t ) 

— directly. 
dT 2 


Here we are faced with the problem of accurately differentiating 
the signals. Only experience and a statistical error analysis 
can tell us which approach is the more accurate, although analy- 
sis by Krause in Reference 12, offers the possibility that sta- 
tistical variation can actually be reduced by proper handling 
of the differentiated signals. 

If G l (Y,7\- l ,T) keeps a constant width as T increases, i.e. , just 
decreases in amplitude, then 


d 2 G, dG, 

— = - 0 where G, = 0 

do 2 da L 
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and we are probably not too wrong in taking 


a 2 G, (Y ,0) / b 2 G, (Y ,0 ,0) / 

-£T / J-T = g e(X> a i>o)/g e (y, 0 ,o) 

Then equation (4.41) becomes 


S(Y,t,0) = M c 





0,a) 



G E ( g>V T > 

g e (y,o,o) 


dT 


(4.42) 

a = 0 

where is an Eulerian time scale. Whether approximations 
like these are valid are only determinable by actual experi- 
mental measurement. 



r 3 2 G L (Y,A 1 , - 0,a) 

IT" 


Homogeneous Turbulence 

Where there is no shear amplification the preceding 
analysis is still valid but the difficulty in taking the fourth 
derivative of the measured covariance becomes even more of a 
problem. It would seem, at this point, that the inaccuracies 
which would result would be difficult if not impossible to 
overcome . 


4. 3. 2. 2 High Speed Jets 

In the Mach emission case 


S(Y,t,0) 


(1 - M cosO)^ 
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(4.29) 



Here, we must measure the curvature of the Eulerian 
covariance or the derivative covariance and integrate along 
the Lagrangian time variable. This can be done as sketched 
in Figure 3. 

The lifetime a* in Lagrangian time must be sufficient- 
ly short that y 2 = A^ sin© = a Q o* sin© is small enough that 
du^/dy^ does not vary appreciably over that distance. If it 
does, then the mean gradient must be left under the integral 
in equation (4.29). 

We have in this case somewhat of a lesser problem in 
taking the second derivative of G, since its curvature is much 
greater than that of G^. Although 3 Gg/dT probably decreases 
as a increases, as long as a* is not too long we can approximate 
equation (4.29) by 


S(Y,t ,0) = - (1 - M cos©) 5 




dT 2 


Y,0,T)1 

J x = 0 


(4.43) 


Then we need only to measure the curvature of the crossed-beam 
autocovariance and the lifetime. Also, we can approximate 
the integrand by 

a 2 G E (Y,0,T)_ 2 ^ 2 G e (Y,A 1 ,T) 

c>T 2 a ° Sa 1 2 

= a 2 G s (Y,A 15 T) (4.44) 

where G is a gradient covariance, 

b 

dp(Y,0,t) d p (Y,0,t + t) 

^y r 

( 12 ) 

and is measurable by crossed-Schlieren techniques. * 
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Figure 3 CURVATURE OF EULER I AN CROSS- COVARIANCE 
REQUIRED FOR MACH EMISSION CASE 
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Then equation (4.43) becomes 


S(Y,t,0) = - (1 - M cosG) 5 aJ o 



G s (Y,0,0) 


(4.45) 


4. 3. 2. 3 General Case 

In the general case we have considered a frame of ref- 
erence (y r , y 2 j y^) with y^ in the direction of the observer 
and left undecided whether the reference frame should be altered 
to move at the speed of sound a Q or the convection speed com- 
ponent a Q M cosQ. The difference between these can be seen more 
clearly by reference to the sketch in Figure 4, showing contours 
of G e (Y, A x , T). 

Assuming that the stress tensor T were relatable to p 

rr r 

(which is not true in general but is at low speeds) , and that 
we can take the time delay derivatives outside of the area in- 
tegral, (which is not justified, section 4.2.3) we are then con- 
cerned with integrals of the form 


r ^ n G(x, 

J dT n 


.^v T) 


dx 


If X_ is taken to move at sonic convection speed, X - a T, 

r * » r i o » 

and we must integrate along lines parallel to = -a Q MT in 
the sketch to satisfy retarded time requirements. 

If X r is taken to move at the particle convection speed 
then we integrate along lines parallel to A^ = a Q MT, to minimize 
spatial gradients. 


From the sketch we can see that integration along 
A r = A r = a Q MT minimizes time scales whereas integration along 
X r = A r - a c T can give very short time scales resulting from 
convection effects. Only at M cos9 =1, as in the sketch in 
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Ai = a o M cosQ T 



& 1 - a o T 


Gg(Y » , T) = const 


)TION OF OBSERVERS REFERENCE 
,ERIAN covariance 
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Figure 5 are both requirements met. At a = 00 , the time delay 
becomes unimportant, (see Figure 5) since changes in T are very 
small. To assess which integration path should be used we need 
a complete knowledge of G( Y, A^., T) or at least an independent 
measurement of gradient fluctuations. 
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for Mach Emission (M cos9 



a Q M cos9 J 


Ajtr a c T) 


for Low Speed Case (a o = 


PATHS FOR M = 1 AND M = 0 


INSTITUTE 


0 M cos9 T 


= 1 ) 


°°) 

CASES 


8 


V. DISCUSSION 

The problems associated with applying crossed-beam tech- 
nology to the measurement of sound source intensities have 
been found to lie in the following general areas: 

1. Measurement of correct covariance 

2. Retarded time effects 

3. Moving frame of reference 

Another area, beyond the scope of this paper involves the 
applicability of point-area correlations and must be proven 
by experiments. Under certain conditions the main problems 
can be circumvented in principle at least. 

Use of the infrared system holds promise of giving measure- 
ments directly relatable to density statistics. This does how- 
ever limit us to sound source integral forms which are expressible 
in terms of density covariances. It has been suggested 
that this representation may result in slowly converging integrals 
and therefore require good measurement accuracy. 

Retarded time and moving axes requirements are made less 
stringent in two special limiting cases, viz., M = 0, and 
M cos9 - 1. In the first of these t the retarded time variable 
can be considered everywhere constant and in the second we 
require only that t = constant in the "plane" of the point- 
plane representation. This latter can be achieved by orienting 
the crossed beams both in the plane of the Mach front. 

The moving frame of reference or Lagrangian covariances 
can be obtained from Eulerian covariances. The main difficulty 
lies in the requirement of obtaining a great many experimental 
measurements to map out fully Gg(A-^,T) and G L (X^,T). The time 
derivatives in the sound source integrals pose another possible 
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restriction in accuracy required of measurement. To obtain 



accurately may require a great many experimental points before 
differentiating. This is alleviated somewhat by the high-shear 
approximation, i.e., determining 



the second derivative rather than the fourth derivative. 

The residual effects of density-gradient fluctuation in 
the convected reference frame may also pose a major problem 
if they are of the same order of magnitude as the "true" time 
derivatives. It is suggested that these be estimated from 
crossed-schlieren results. 

An attempt to obtain a unified form of the sound source 
integral applicable to crossed-beam measurements at all 
speeds has been unsuccessful here. Opposing requirements on 
retarded time and space-derivative minimization are a major 
factor in this conclusion. Also, representation of the stress 
tensor in terms of a parameter (like density) which can be 
measured by crossed-beam methods appears to be impossible at 
other than M = 0 and M cosG = 1. However it is believed 
that measurements using the M = 0 and M cos9 = 1 approximation 
should cover nearly the complete speed range of interest. 

Most serious of all assumptions made for low speed studies, 
is the assumption that the jet is isothermal. The range of 
validity of this approximation must be determined. Also 
further analysis on entropy fluctuation measurement must be 
made. 
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